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1.  Introduction 


The  free-flight  motion  of  a  slowly  spinning  elastic  missile  has  been  studied  by  a  number  of 
authors  (1-6).  Although  the  theory  of  Murphy  and  Mermagen  (6)  is  valid  for  large  spin  rates,  no 
Magnus  force  is  included  in  the  aerodynamic  force  distribution.  In  this  report,  the  force 
distribution  is  extended  to  include  a  Magnus  force  distribution  function.  The  finite  element 
method  (FEM)  (6)  is  then  used  to  calculate  the  first  five  positive  frequencies  and  the  first  three 
negative  frequencies  of  the  motion  of  a  10-cal.  spinning  cone  cylinder.  The  elastic  frequencies 
are  shown  to  be  strongly  affected  by  the  high  rate  of  spin  required  to  stabilize  this  projectile.  For 
certain  values  of  the  elasticity,  it  is  shown  that  spin  resonance  can  occur  with  the  lower  elastic 
positive  frequencies,  and  maximum  strain  can  exceed  the  plastic  limit. 

2.  Coordinate  System 

The  elastic  missile  is  assumed  to  consist  of  a  very  heavy  elastic  circular  rod  of  fineness  ratio,  L, 
and  mass,  in,  embedded  in  a  very  light  symmetric  aerodynamic  structure  that  may  be  longer  than 
the  rod.  The  rod’s  axial  moment  of  inertia  is  Ix ,  and  its  transverse  moment  of  inertia  about  its 
center  is  Il0 .  The  rod’s  diameter  can  vary  over  its  length,  and  its  maximum  diameter  will  be 
denoted  by  d.  All  distances  will  be  expressed  as  multiples  of  the  rod  diameter,  and  its  length  is 
Ld.  A  nose  windshield  of  length,  x2id ,  may  be  attached  to  the  forward  end  of  the  rod,  and  fins 
or  a  boat  tail  may  extend  beyond  the  end  of  the  rod  at  a  distance,  x01d.  Thus,  the  rod  is  located 
between  x1  =  -L/2  and  x2  =  L/2  ,  while  the  aerodynamic  structure  extends  from  x0  =  x,  -x01  to 
x3  —  x2+  x23 . 

An  earth- fixed  coordinate  system  will  be  used  with  the  X(-axis  oriented  along  the  initial  direction 
of  the  missile’s  velocity  vector.  The  Ze- axis  is  downward  pointing  and  the  Yc.-ax  i  s  detennined  by 
the  right  hand  rule.  A  nonrotating  coordinate  system,  XYZ  is  then  defined  with  origin  always  at 
the  center  of  the  rod  and  the  X-axis  tangent  there.  The  X-axis  pitches  through  the  angle,  6 ,  and 
yaws  through  the  angle,  y/,  with  respect  to  the  Xe-axis.  Body-fixed  coordinates,  XYhZh ,  are  now 
defined  for  which  the  Y/-Z/,  axes  rotate  with  the  missile. 

We  will  conceptually  slice  the  missile  into  a  large  number  of  thin  disks  perpendicular  to  the 
X-axis  with  thickness,  dx.  When  the  rod  flexes,  the  disks  shift  laterally  perpendicular  to  the 
X-axis  and  cant  to  be  perpendicular  to  the  centerline  of  the  disks.  This  canting  action  neglects 
the  shear  deformation  of  the  rod,  and  this  constraint  is  called  the  Bernoulli  assumption  (7).  The 
lateral  displacement  of  a  disk  has  body-fixed  coordinates  5 hv ,  5 bz ,  and  the  disk  is  canted  at  angles 

r  r 

±  by  ’ 1  bz  • 


1 


(1) 


r, 


M 


by  . 


by 


dx 


r, 
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dx 


It  is  important  to  note  that  at  the  central  disk 


4  (o.0  =  4  (°.0  =  0,  (o.O = r„  (o,/) = o . 


(2) 


The  earth-fixed  coordinates  of  the  central  disk  are  (xe,  ye ,  ze ),  and  the  earth- fixed  coordinates  of 
the  other  disks  are  computed  in  terms  of  the  central  disk  earth-fixed  coordinates,  their  body-fixed 
displacements,  and  the  Euler  angles  9,yj,(j)\ 


and 


Xde=Xe+X 


1-1 


(v2  +02)/2  ]-wRe{(Sby  +  iSb:y}  +  9Im{(Sby  +  id^e*} 


yde  =ye+W  +  Re[[dhy+  iShz )  e' ’*  } , 
=ze-xd  +  Im { (shy  +  iShz )e¥J  . 


(3) 

(4) 

(5) 


Murphy  and  Mennagen  (5,  6)  used  the  nonspinning  elastic  coordinate  system  with  XYZ  axes. 
(0  =  0) .  The  lateral  displacements  of  a  disk  in  this  elastic  coordinate  system  are  shown  in 
figures  1  and  2  and  can  be  computed  from  body-fixed  quantities. 
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Figure  2.  X-Y  coordinates  of  the  cross-sectional  disk. 


SE  =  SEy  +  iSEz  =  (Sby  +  *<*fc  )  ^  ,  (6) 

and 

r  =  r>  +  iT,  =  (rby  +  irky.  (7) 

In  Murphy  and  Mennagen  (5),  the  partial  differential  equation  (PDE)  for  the  missile’s  flexing 
motion  is  derived  by  use  of  Newtonian  mechanics  and  the  assumption  that  each  disk  can  be 
assumed  to  be  a  point  mass.  Thus,  the  canting  of  the  disks  is  neglected,  and  their  axial  and 
transverse  moments  of  inertia  are  set  equal  to  zero.  The  complete  Lagrangian  is  obtained 
including  the  canting  of  the  disks  required  by  the  Bernoulli  assumption  and  is  used  to  obtain  a 
more  accurate  PDE  in  appendix  A.  For  the  small  spin  considered,  this  modification  has  a  very 
small  effect  on  the  eigenfrequencies.  For  a  spin-stabilized  projectile,  the  required  spin  can  be 
30-100  times  the  aerodynamic  frequencies  and  can  have  a  large  effect  on  the  eigenfrequencies. 


3.  Aerodynamic  Force 


The  ye  and  z(,  components  of  the  central  disk  velocity  can  be  approximated  by  linear  relations  in 
angles  of  pitch  and  yaw  with  respect  to  inertia  axes  (6, if/)  and  angles  of  attack  and  sideslip  with 
respect  to  the  velocity  vector  (a,/?)  (figures  1  and  2). 
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and 


ye=(vld){P+v)’ 


(8) 


ze={V/d)(a-e).  (9) 

Equations  8  and  9  can  be  written  as  a  single  complex  equation: 

ye+iZe={v/d){<li+<lie)’  0°) 

where 

qt  =  fi  +  ia  (11) 


and 


qu  =V~iO. 


(12) 


In  Murphy  and  Mermagen  (5),  the  linear  aerodynamic  force  loading  is  expressed  in  terms  of 
three  force  distribution  functions,  cD  [x),cfI  (x),c/2  (x)  and  the  base  pressure  coefficient,  CDb  , 
plus  a  body-fixed  force  associated  with  possible  bent  fins.  The  complex  angular  velocity  of  the 
central  disk  is  q2  =  qle .  Because  the  lateral  motion  of  the  missile  is  quite  small,  qx  =  -q2 ,  and 
the  aerodynamic  damping  force  terms  in  the  aerodynamic  loading  on  the  aerodynamic  structure 
can  be  combined.  For  a  rapidly  spinning  projectile,  we  will  introduce  a  Magnus  force 
distribution  function  c  fM  (x) .  This  aerodynamic  loading  in  nonrotating  elastic  coordinates  is 

dF 

-jL  =  -gfD{x)-  O3) 

dx 


dFv  .  dF 

-  +  l  - 


dx  dx 


Si 


;fl(x)  +  i(<f>d/v)cJM(x)  q^T  +  (SE  - xq^d/V) 


+  c 


f  2 


(x)(2qt-t)(dlV) 


(14) 


Fxhp  =  -SiCt 


Dbp  ' 


(15) 


The  total  aerodynamic  force  acting  on  the  aerodynamic  structure  is  given  by  the  integrals  of 
equations  13  and  14  and  by  adding  the  base  drag  of  equation  15  to  the  axial  force: 


and 


F  =  Si 


[cx+i(<j)d  lv)cXM^qx+c2(qxd /V) 

~ J\{t)-i(<f>d/v)jm(t)-j2(t)(d/V ) 


(16) 


(17) 


where  various  functions  are  defined  in  appendix  B. 
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Similarly,  the  transverse  aerodynamic  moment  about  the  rod’s  center  can  be  computed  from  the 
transverse  aerodynamic  force  and  a  small  axial  force  contribution: 


M  =  My  +  iMz 


(<c3+i^d/v)c3M)ql+c4(qld/V) 

-J3(t)-i{^d/v)J3M(t)-j4(t)(d/v)-J5(t) 


(18) 


The  primary  components  of  drag  are  head  drag  and  base  pressure  drag.  The  third  component  is 
skin  friction  drag  that  is  -15%  of  the  total  drag  and  will  be  neglected  in  this  report  to  simplify 
the  FEM  calculations. 

The  actual  angular  motion  will  be  described  by  aerodynamic  moment  about  the  center  of  mass. 
The  static  and  Magnus  moment  coefficients,  for  example,  are  as  follows: 

C Ma  _  C3  _  XcC\  ’  (19) 

and 

C Mpa  ~  C3 M  ~  XcC\M  ’  (20) 

where 

X2 

xc  =  (1/Z)  jxpjbx  .  (21) 

*i 

The  center  of  mass  of  the  projectile  will  be  at  the  center  of  the  rod  if  the  axial  mass  density  along 
the  rod,  px ,  is  constant. 


4.  Frequencies 


For  simplicity,  the  spin  will  be  required  to  be  always  positive.  The  angular  motion  of  a  statically 
stable  rigid  missile  ( CMa  <  0 )  can  be  described  by  the  sum  of  two  complex  exponentials: 


Kl0e 


!  *  ,  K  S*2+d> 2, 


(22) 


where 


and 


1 1  £  ]  d(--Ma 


(24) 
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These  two  frequencies  are  opposite  in  sign  and  <j>XR  >  <f>2R  .  The  motion  of  a  statically  unstable 
missile  ( CMa  >0)  will  have  the  same  expression  if  it  is  gyroscopically  stable  (sg  >  1 ).  (j>lR  and 
<j)2R  are  both  positive,  but  c/)]R  is  still  larger  than  <f>2K . 

For  a  nonspinning  statically  stable  missile, 

</>lR  if>2R  ~  coR  ,  (25) 

where 


The  frequencies  for  a  statically  unstable  missile  with  gyroscopic  stability  factor  equal  to  1  have  a 
similar  fonn. 

fi\R  ~  02R  ~  R  '  (27) 

The  elastic  motion  of  a  homogeneous  circular  rod  with  constant  diameter  is  detennined  by  its 
fineness  ratio  L.  and  an  elastic  parameter  col . 

col  =  E0I0L/md 3 .  (28) 

E0  is  Young’s  modulus  at  the  center  of  the  rod  and  I0  is  the  area  moment  of  inertia  at  the  center 
of  the  rod.  The  standard  analysis  for  a  nonspinning  free-free  beam  gives  the  following  relation 
for  infinity  of  elastic  frequencies: 

(29) 

where 

fK  =  4.730,  7.853, 10.996, 14.137...  .  (30) 

For  the  odd  number  modes,  the  rod  is  symmetric  U-shaped,  whereas  for  the  even  number  modes, 
the  rod  is  antisymmetric  S-shaped. 

In  Murphy  and  Mennagen  ( 6 ),  the  parameter  a  =  coxlcoR  was  used  as  a  measure  of  the  elasticity 
of  a  finned  missile.  We  will  continue  to  use  this  parameter  to  describe  the  elasticity  of  a  spin- 
stabilized  projectile.  For  a  slowing  spinning  finned  projectile,  it  has  been  shown  that  the 
aerodynamic  frequencies  are  affected  by  the  elasticity  when  a  <  20  .  For  a  spin-stabilized 
projectile,  the  first  elastic  frequency  should  be  compared  with  the  spin  and  not  the  aerodynamic 
frequency.  Thus,  we  will  see  that  the  elastic  frequencies  are  affected  by  the  spin  when  cr  <  200  . 

The  frequencies  present  in  the  motion  of  a  nonspinning  elastic  projectile  would  form  an  infinite 
sequence  where  the  first  two  frequencies  would  be  related  to  cj)XR  and  (j)2R  ,  while  the  later  ones 
would  evolve  from  ±coK  ,  i.e.,  [<p2K+]  =  coK;<j)2K+2=  -coK  for a>  2 00) .  The  odd-numbered  modes 
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rotate  in  the  direction  of  the  spin,  have  positive  frequencies,  and  are  called  positive  inodes,  while 
the  even-numbered  modes  have  negative  frequencies  and  are  called  negative  modes.  Since  the 
modes  for  a  spinning  missile  essentially  bifurcate,  j  =  3  is  a  symmetric  mode  rotating  in  the 
direction  of  spin,  j  =  4  is  a  symmetric  mode  rotating  in  an  opposite  direction  to  the  spin,  j  =  5  is 
an  antisymmetric  mode  rotating  in  the  direction  of  spin,  j  =  6  is  an  antisymmetric  mode  rotating 
against  the  direction  of  spin,  etc. 


5.  FEM 


The  rod  is  assumed  to  be  represented  by  the  sum  of  an  inelastic  bent  component  rotating  with  the 
missile  and  an  elastic  deformation. 

dE{x,t)  =  8EB(x)e* +dE{x,t\,  x,<x<x2,  (31) 

and 

db(x,t)  =  dEB(x)  +  db(x,t)-,  x,<x<x2,  (32) 

where 

,  .  ddFR( 0) 

<W°)  =  -fp  =  °-  03) 

Because  the  aerodynamic  nose  structure  is  rigidly  attached  to  the  rod, 

§E  (x,t)  =  §E  (x2,t)  +  (x-x2)  r(x2,t)  x2<x<x2.  (34) 


The  motion  of  the  elastic  component  of  the  rod  is  controlled  by  the  elasticity  of  the  rod  and  the 
aerodynamic  force  acting  on  it. 

FEM  is  a  very  powerful  method  for  calculating  the  time  history  of  the  elastic  flexing  motion. 

The  rod  is  divided  into  nj  elements  of  length  Le  =  Ljrij .  We  will  consider  only  an  odd  number 
of  elements  with  the  center  of  the  central  element  satisfying  equation  2. 

The  shape  of  the  j-th  element  is  given  by  a  linear  combination  of  third-order  Hermitian 
polynomials  (6). 

db(x’t)  =  y£jqbp(t)Np(z),  (35) 

1 

where 

x  =  Le(z  +  Zj ),  (36) 

Zj=Xl/Le+J~l  ’  (37) 

and 

0  <  z  <  1 .  (38) 
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The  coefficients  of  the  polynomials  are  complex  functions  of  time  and  are  called  connectors. 
The  first  two  connectors  are  the  deflection  and  slope  of  the  left  end  of  the  element  and  the  third 
and  fourth  are  the  deflection  and  slope  of  the  right  end.  To  ensure  continuity  in  deflection  and 
slope  at  junction  points,  the  corresponding  pairs  of  connectors  are  equal.  For  n  j  elements,  there 
are  2 «  independent  complex  connectors,  qn .  It  is  convenient  to  let  the  index  for  the  connectors 
run  from  3  to  n,  =2ni+2. 

In  Murphy  and  Mermagen  (6),  n,  complex  second-order  differential  equations  are  derived  for 
the  nt  complex  variables  qn . 


+(5'»m  +  it/>Smn)^n  +  (Tmn  +  ^Tmn  )  Cl, 

n  =  1 


(39) 


where 

=  (yhny  +  / 


(40) 


and 


tm=(md2)tmD+(gld)tmA 


(41) 


The  In,  coefficients  in  equation  41  are  defined  for  no  Magnus  force  in  Murphy  and 
Mermagen  (6).  The  inclusion  of  the  Magnus  force  distribution  function,  cfM  (x) ,  modifies  some 
of  the  T*m  and  tmA  coefficients  in  appendix  C. 


6.  Cone  Cylinder  Frequencies 


Transient  frequencies  and  damping  rates  for  a  9-cal.  cylindrical  rod  with  a  1-cal.  conical  nose 
can  be  obtained  from  the  homogeneous  part  of  equations  41,  ( tm  =  0  ).  The  necessary  parameters 
are  given  in  appendix  D.  The  first  eight  frequencies  and  damping  rates  for  a  =  20  were 
obtained  from  3-,  5-,  and  7-element  codes  and  are  compared  with  results  of  the  PDE  method  of 
Murphy  and  Mennagen  (5)  (table  1).  The  frequencies  that  differ  from  PDE  results  by  >5%  are 
marked  by  a  “x.”  The  3-element  code  gives  good  results  for  the  first  four  frequencies  and  the 
5 -element  code  is  good  for  the  sixth,  seventh,  and  eighth  frequencies.  The  fifth  frequency, 
however,  requires  more  than  seven  elements.  In  figure  3,  the  rod  shapes  for  the  fifth  mode  are 
plotted  for  the  FEM  code  and  the  PDE  value.  We  see  that  the  antisymmetric  shape  specified  by 
the  PDE  value  probably  requires  at  least  13  segments  to  describe  it. 
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Shape  of  the  beam  at  sigma  =  20,  p  =  7000 


Solid  Curve  =  PDE 

- (dash  dash)  Curve  =  3  element 

(dash  dot)  Curve  =  5  element 
. (dot  dot)  Curve  =  7  element 


Figure  3.  Shape  of  the  beam  for  j  =  3,  5,  7  element  and  PDE  (real  part  of  the  shape). 

The  two  aerodynamic  frequencies  do  not  differ  from  their  rigid  values.  The  elastic  frequencies, 
however,  are  quite  different  from  nonspinning  values  of  ±eoK  .  The  first  two  positive  frequencies 
are  150%  and  130%  greater  than  these  values. 

Figure  4  shows  the  positive  and  negative  first  elastic  frequencies  divided  by  the  zero-spin  value 
given  by  equation  30.  When  the  first  elastic  frequency  is  much  greater  than  the  spin  ( cr  «  200  ), 
these  ratios  are  near  unity.  Figure  5  shows  similar  results  for  positive  and  negative  second 
elastic  frequencies. 


7.  Bent  Projectile  Resonances 


In  Murphy  and  Mermagen  (5,  6),  a  bent  rod  was  described  by  a  pair  of  quartic  curves: 

-2  '  J  ~4  -L/2  <  x  <  0 

2 


8 EB  — 1/|  |X  +  of,  jX 


=  dl2x  +  d21x 


0  <  x  <  L/2. 


(42) 
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hi5dot/w2,  -phi6dot/w2 


Figure  5.  (j)5  j  (02,~  (j)(}  jco2  vs.  (J. 
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These  curves  can  then  be  used  to  calculate  the  inhomogeneous  coefficients  in  equation  42.  For  a 
slowly  spinning  finned  missile,  the  aerodynamic  coefficients,  tmA ,  dominate  while  for  a  rapidly 
spinning  projectile  the  dynamic  coefficients,  tmD ,  dominate.  These  inhomogeneous  tenns  induce 
a  modal  response  at  the  spin  frequency  that  can  have  a  large  amplitude  when  the  spin  is  equal  to 
one  of  the  transient  frequencies.  The  values  of  d(>.  given  in  appendix  D  specify  the  very  small 
deflection  of  the  rod  forward  tip  of  0.0015  in. 

For  constant  spin,  special  solutions  of  equation  39  have  been  calculated  in  Murphy  and 
Mermagen  (6).  These  trim  solutions  have  the  form 

(43) 

Using  5  elements,  the  12  sn ’s  can  be  computed  for  fixed  values  of  spin  and  cr .  The  complex 
location  of  the  forward  end  of  the  rod  is  specified  by  vn . 

In  figure  6,  the  amplitude  of  the  forced  motion  of  the  forward  end  of  the  rod  is  plotted  vs.  cr  for 
two  values  of  spin.  The  resonances  at  cr  =  1 19.3, 126.3  occur  when  the  positive  first  elastic 
frequency  is  equal  to  a  spin  value  while  the  resonances  at  cr  =  21.9,  23.2  occur  when  the  positive 
second  elastic  frequency  is  equal  to  a  spin  value.  The  amplitude  of  the  flexing  motion  of  the 
forward  end  of  the  rod  is  70%  of  its  diameter. 
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In  Murphy  and  Mermagen  (6),  a  nonlinear  spin  equation  is  derived: 

1 A  +  Re{md 1  (Qn  -  iq2J6  +  iq2J, s )  -  gxd  QAJ  =  gxd(Cx  )w ,  (44) 

where  the  cubic  terms  IX,QD,JS,  and  QA  are  defined. 

For  a  cone  cylinder,  the  linear  roll  moment  coefficient  has  a  very  simple  form, 

(c4w=c,,(^M-  <«) 

Equations  39  and  44  for  five  elements  can  be  integrated  to  show  the  occurrence  of  resonance 
with  the  positive  first  elastic  frequency  (er  =  122.8).  All  initial  conditions  are  made  zero  except 
for  (j){]  =  7000  rad/s.  Figures  7-9  show  the  time  variation  of  spin,  angle  of  attack  magnitude,  and 
rod  forward  tip  motion  magnitude,  <f),  | qxb  | ,  \qxxb  | .  Resonance  is  clearly  shown  at  t  =  0.25  s  when 
spin  is  near  7000  rad/s.  According  to  figure  9,  the  maximum  amplitude  of  the  forward  rod 
motion  is  30%  of  the  rod  diameter.  Thus,  the  motion  amplification  due  to  spin  going  through 
resonance  is  less  than  half  its  resonance  value.  For  rigid  finned  projectiles,  relations  between 
pitching  motion  amplitude  due  to  spin  varying  through  resonance  and  its  resonance  value  are 
given  in  Murphy  ( 8 ). 


phidot  versus  time 


Figure  7.  (j)  vs.  time  for  <j=  122.8,  =7100. 
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mag(qllb) 


Figure  8.  mag[qlb )  vs.  time  for  <r=  122.8,  (f)Q  =  7100. 


Figure  9.  mag(qllh )  vs.  time  for  cr  =  122.8,  (j)0  =  7100. 
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For  the  symmetric  waveform  of  the  first  elastic  mode  for  a  nonspinning  rod  with  no  force 
loading,  the  maximum  strain  occurs  at  the  midpoint. 


(1/2) 


>4(  0) 


A 


(46) 


Figure  10  shows  the  variation  of  maximum  strain  with  time.  For  most  metals,  yield  occurs  for 
strain  >0.0015.  Thus,  figure  10  shows  yield  at  resonance. 


mag(eM)  versus  time 


Figure  10.  mag  (eM  )  vs.  time  for  a  =  122.8,  =  7100. 


8.  Summary 


The  previously  derived  FEM  theory,  which  was  applied  to  slowly  spinning  finned  projectiles, 
has  been  extended  to  rapidly  spinning  spin-stabilized  projectiles. 

Positive  and  negative  elastic  frequencies  for  a  spin-stabilized  projectile  have  been  calculated,  and 
their  magnitudes  have  been  shown  to  be  significantly  different  for  a  <  200  . 

Resonance  with  the  first  elastic  mode  has  been  demonstrated,  and  very  small  rigid  asymmetries 
have  been  shown  to  cause  yield  at  resonance. 
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Appendix  A.  Improved  Partial  Differential  Equation 


The  roll  moment  of  inertia  of  a  circular  disk  is  2ad  (md2  ]p2(x  )dx  and  its  transverse  moment  of 
inertia  is  ad(md2  ]p2(x)dx  where  ad  is  ( 1 6  L )  'and  p2  (x)  describes  the  variation  of  moments  of 
inertia  along  rod.  In  the  derivation  of  the  partial  differential  equation  (PDE)  in  Murphy  and 
Mermagen,1  the  disks  were  assumed  to  be  point  masses,  i.e.,  ad  =  0 .  In  Murphy  and 
Mermagen,2  the  kinetic  energy  associated  with  ad  for  a  disk  was  shown  to  be  Taddx  ,  where 

Tad=(admd2/2)[f-;+f2  +  2Re{(i)(2if-(i)r)f  +  i(f-2ici)r)Q}]p2,  (A-l) 

where  Q  =  iqle . 

According  to  Geradin  and  Rixen,3  two  terms  in  Tad  appear  in  the  PDE  for  a  flexing  projectile: 


d  dT 


dt  ar. 


=  (md2ad )  [f  y  +  4<j)f  z  +  2^Ty  +  r  -  2<j>q]  , 


(A-2) 


d  f  8T 


dt  dr. 


=  ( md2ad )  [f  z  -  4<j)ry  +  2<j)2rz  -  q  -  2<j>r] 


(A-3) 


The  contribution  of  Td  to  the  improved  PDE  is  obtained  by  multiplying  the  second  term  by  i, 

adding  it  to  the  first  tenn  and  differentiating  the  result  with  respect  to  x.  The  improved  version 
of  equation  52  in  Murphy  and  Mermagen1  is 


d~SE  2d4SE  2f 
dt2  +(°0  dx4  +C0°kCd4 


Si  L(cn  +i(pd/V)cfM)r  +  (c/2f  -cnSE)(d/V ) - 


adL—(  f  -  4/^r  +  2^2r )  +  E£  +  E2£  +  ixQ  -  N. 


(A-4) 


1  Murphy,  C.  H.;  Mermagen,  W.  H.  Flight  Motion  of  a  Continuously  Elastic  Finned  Missile.  Journal  of  Guidance,  Control 
and  Dynamics  January-February  2003,  26,  89-98. 

9 

“  Murphy,  C.  H.;  Mermagen,  W.  H.  Spin-Yaw  Lockin  of  an  Elastic  Finned  Projectile ;  ARL-TR-3217;  U.S.  Army  Research 
Laboratory:  Aberdeen  Proving  Ground,  MD,  August  2004. 

3 

Geradin,  M.;  Rixen,  D.  Mechanical  Vibrations:  Theory  and  Applications  to  Structural  Dynamics',  John  Wiley:  Chichester, 
1994,  pp  172-174. 


17 


The  first  boundary  conditions  at  Xi  and  X2  are  modified  by  adding  the  sum  of  the  first  term  and  i 
times  the  second  tenn  evaluated  at  x  1  and  X2.  Equations  53  and  55  of  Murphy  and  Mermagen1 
become 


d3SE(xl,t ) 


dx 


+  kcd3(xl,t)  +  g4cD(xl)S(xl,t) 


+  adL 


iQ~{y  (x1,?)  +  2^2r(x1,t))-2^(ig-2f  ( 


xx,t 


CO  n  =- 


gJi  (A-5) 


and 


83SE(x2,t) 


dx 3 


+  kcdi(x2,t)  +  g4cD(x2)S(x2,t) 


+  adL 


iQ-(j'(x2,t)  +  2^2T{x2,t))-2i^(iQ-2t{x1,t))  co02  =  g4f2.  (A-6) 


The  trim  solution  equations  78,  79,  81,  86,  and  87  of  Murphy  and  Mermagen1  are  replaced  by 

dcn 


d  dET  d  8 ET  dS ET 

dx 4  6  dx2  4  dx 

=  E^t-N*t-E*bf  +  Eb, 


E5Set  g4dTl 


dx 


(A-7) 


and 


where 


+  g4c0  (x,  )<Jr  (*, )  +  E„  [(4/5)  -  rr  (*, )] 

—  —  LNTco0~  J  —  ~g4f\T, 

7AM  +  gtcD(x2)Sr  (x2 )  +  E„  [(4/5)  -  rr  (x2 )] 
(^/F)[g4Cia4  -MVv]  =  g4/2r, 


dx 


.4 


dx 


2  -^4 


U'l 


dx 


+  E5w 1 


^3<=’7’  Ay  EBF  +  Tig  +  g4dcT 


dc 


D 


dx 


(A-8) 


(A-9) 


(A- 10) 


dwm 

dx 


+  E<wm=  0; 


m  =  2,  3  , 


(A- 11) 


a4 

a6 


&4 


C/1  + 


C/2  +  C/M  )  CD 


=  5 adL<f>2co(  2 


0  ’ 


(A- 12) 
(A- 13) 


18 


and 


Eb=8a 


(cfl  +i(*d/v)c/M)rB+(cf2 Tb  - cflSEB)(i</>d/v)-d^ ^ 


(<j>/(0°)2  SB+E6^T 


(A-14) 


The  transient  solution  equations  92,  94-96,  and  98  of  Murphy  and  Mermagen1  are  replaced  by 


and 


where 


and 


dVk  tj  dVk  c  dv|/k 


dcr 


b,  -  E.  -  E4k  ^  +  E5kv|/k  -  g4ykc  —  =  N: 


dx4 

d4w. 


dx4 
d4w„ 


-E 


c  dx2 
d2w. 


dx 

dw,( 


dx 


"3k  x,k’ 


dCr 


dx 


6k  ,  2  -E4k-rL  +  E5kW4=E3k-Nk+g4Vkc 

dx  dx  dx 

d2w  dw 

-  E„. - -  •  H  „ — — +  E5kwm  =  0  i  m  =  5,6, 


"6k 


dx" 


"4k 


dx 


-yy-  +  §4E- nCD  ( wk  -  Vkc )  +  V 


dx 


Elk  E6k 


d\j/[xx ) 


dx 


a,  (At  - 2i<f,)(d/V)[g4Cu  - LN^bl'  =  -gtE;lfm 


A  +  sXXd  ( wk  -  Vu, )  +  K 


di//(x2) 

dx 


+  ad  ( A k  -2i<j>)(d/V)[g4Cla -LN.coy^b'1  =  gAE~lf2k, 


Ee k  =  adL  (4  -  4i<f>Ak  +  2<t>2 )  ®o2 > 

Elk  =  ad^Ak  ( A  _  )  0J(>  , 


63  =  1  +  2ko)x  1  ( Ak  -  itff'j . 


(A- 15) 
(A- 16) 
(A- 17) 

(A- 18) 


(A- 19) 

(A-20) 

(A-21) 

(A-22) 


For  simplicity,  the  small  Magnus  contribution  to  the  boundary  conditions  in  this  appendix  have 
been  neglected. 
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Appendix  B.  Integrals 


B.l  Aerodynamic  Coefficents 

*3 

c)  =  f  cfldx  =  CNa 

x0 

_ © 

a 

ii 

•3 

o 

II 

CO 

x3 

c\  m  ~  j*  cjMdx  =  CNpa 

Xq 

x3 

C3M  =  J  CfM XC^X  ~  [^M/?a]0 

x0 

x3 

C2=\CfldX=\CNq+Cm\ 

Xq 

x3 

c4  =  J  cf3xdx  =  [CM?  +  CMd  ]Q 

x0 

cf3  =  2c  f2  —xCfX 

[  ]o  ^  Xc  =  0 

B.2  Functions  of  Time 

5c(0  =  (1/Z)J5£A^ 

*1 

oS 

II 

o 

•-s 

X 

8* 

> 

II 

o 

■-s 

8* 

-'■3 

Jm  (0  =  j  c  fM  r  xdx 

Xq 

ll 

o 

* 

8- 

x3 

J 4  (V)  =  J  [cfir —  Cffi^xdx 

Xq 

1 

(N 

II 

"-T 

J5  0 t)  =  J  CD  (< -  8C )  dx  +  [c S£  (. Xy )  -  8C ]  CDbp 

x0  *0 
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Appendix  C.  Magnus  Terms 


x3 

x3 

cim  —  CNpa  —  J  cfM  (x)  dx 

xo 

C3M  —  CMpa  —  J  cfM  (x)xdx 

xo 

x3 

^imb  ~  J  cairBdx 

xo 

x3 

^3MB  _  J  CfM^BXdX 

xo 

^qj  =lcfM(X)N;(z)dz 

0 

x  =  Le(zj+z) 

K  =  j[c™(x)NUz)]xdz 

0 

Zj  =(Xl/Le)  +  j-1 

1 

hpji  =Lejc)M(x)Np(z)dz 

0 

p,q  =  l,2,3,4 

1 

^Bpj  —  Lgjc^FgNpdz 

0 

^  qc  9  ^  qa  ’  ^  qa  9  ^  qc  9  ^  pci 9  ^  pal  9  ^  Bpc 

,h„  are  calculated  from  h ,  h  2. 

opa  qj  '  qj 

same  manner  as  /  ,  /  were  computed  from  / . . 
hin,h2n,hml ,hBm  are  then  computed  in  the  same  way  as  fm  . 


^  Pqj  =  |cfM  (x)Nq  (z)Np  (z)dz  x  =  Le(zj  +  z) 

0 

Zj  =(Xl/Le)  +  j-1 

p,q  =  l, 2,3,4 


hm  is  computed  from  hpqj  in  the  same  manner  as  fmn  was  computed  from  /  . . 
m,  n  >  3  . 
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T;,=(g,dj/V)c,„ 

T;,=-(g,dVv)c!M 

t;2 = -ix 

T,:=-(g,d2/v)h,„ 

T;n=-(g,d!/v)h,n 
Tli  =-(g]d2/v)hml 

^  =  2md2adb2m 

T1  =-(gld2/V)hmn-(mdVL)(2®o/®l)kCran  111,  II  >  3 

tlA  =  JlB+i(<M/V)(J2B+JlMB) 

t2A  =  J3B+i(<i)d/V)(J4B+J3MB) 

U  =  fB,n  +  faBm  +  1  («j>d/V)  (gBm  +  SaBm  +  hBm  )  ™  >  3 


The  Magnus  force  on  the  aerodynamic  extension  has  a  small  effect  on  the  flexing  motion  and  has 
been  neglected  in  the  previous  relations. 
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Appendix  D.  Projectile  Parameters 


P\  P2  P2  1 


x  =  0 


L  =  9 


V=  3000  ft/s 


d  =  0.35  ft 


p=  0.002  slugs/fr 


m  =  1 .60  slug 


*01  =  0  *23  =  1 


Ix=  0.0245  slug-fr 


It  =  1.335  slug- ft" 
Cf J  =4(5.5  -  X ) 

_  2(x-5.5) 


4.5  <  x  <  5.5 
-4.5  <  x  <  4.5 


cf2  =  2(5.5  -x) 

=  2  +  0.5(l-e2M‘ 


4.5  <  x  <  5.5 
-4.5  <  x  <  4.5 


C./M  =  0 


=  -0.07 


=  (0.30X5.5  -x) 


2  <  x  <5.5 
-  4.5  <  x<2 

4.5  <  x  <  5.5 
x  <  4.5 


cDbp  =0.14 


C,„  =  -0.025 


dl2=  2  x  10 


d22  =-1.25  x  10 


</>  =  7000  rad/s 
co.  „  =  102.8  rad/s 


co^  =  5 1.3  rad/s 

C  =2  5 
CNq  +  CNa=  34 
C*.=-0.5 


d\\  ^^21  ^ 

co2i?  =  25.6rad/s 

sg  =1-6 

cMa  =  n.i 

CMq+CMa=-  50 

cw.=u 
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a  =  20  <j)/coR  =  136  coRl/coR=  1.99 
co1/coR=  20.0  co2/coR=  55.1 

Table  D-l.  Transient  frequencies  and  damping  rates. 


cor2/cor=.50 
co3/coR  =108 


Code 

k 

A 

\  /  Mr 

3  element 

i 

2.003 

-0.0074 

5  element 

i 

2.002 

-0.0074 

7  element 

i 

2.003 

-0.0074 

PDE 

i 

1.993 

-0.0075 

3  element 

2 

0.502 

-0.0263 

5  element 

2 

0.502 

-0.0263 

7  element 

2 

0.502 

-0.0263 

PDE 

2 

0.505 

-0.0265 

3  element 

3 

51.25 

-0.0113 

5  element 

3 

50.91 

-0.0168 

7  element 

3 

50.85 

-0.0167 

PDE 

3 

50.71 

-0.0112 

3  element 

4 

-32.85 

-0.0213 

5  element 

4 

-32.76 

-0.0139 

7  element 

4 

-32.74 

-0.0138 

PDE 

4 

-32.67 

-0.0139 

3  element 

5 

138. 61x 

-0.0186 

5  element 

5 

130. 02x 

-0.0274 

7  element 

5 

125. 93x 

-0.0256 

PDE 

5 

118.94 

-0.0237 

3  element 

6 

-70.48 

-0.0305 

5  element 

6 

-69.20 

-0.0161 

7  element 

6 

-68.57 

-0.0159 

PDE 

6 

-67.29 

-0.0158 

3  element 

7 

172. 74x 

-0.0271 

5  element 

7 

161.55 

-0.0255 

7  element 

7 

160.45 

-0.0245 

PDE 

7 

159.76 

-0.0326 

3  element 

8 

-01.99x 

-0.0164 

5  element 

8 

-93.05 

-0.0150 

7  element 

8 

-92.62 

-0.0148 

PDE 

8 

-92.34 

-0.0144 

Note:  PDE  =  partial  differential  equation. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


cfj  M 
d 

E(x) 

E0 

F 

ad 

8i 


I(x) 

h 

lx 
1, 0 
L 

4 

m 

P  =  j> 

nj 

n, 


di 

die 


<h 

d„ 

dbn 


aerodynamic  force  distribution  functions 

maximum  rod  diameter 

Young’s  modulus 

Young’s  modulus  at  rod  center 

F  +  iF,  complex  transverse  aerodynamic  force 

(161)-' 
pV2S/ 2 

(df  jj  ydydz  =  (df  JJ  z2dydz  ,  area  moment  of  rod 

area  moment  at  rod  center 

axial  moment  of  inertia  of  projectile 

transverse  moment  of  inertia  of  projectile  about  rod  center 

rod  length/rod  maximum  diameter 

L/rij ,  dimensionless  length  of  element 

projectile  mass 
spin 

number  of  rod  elements 

2  rij  +2 

P  +  ia  ,  complex  angle  of  attack  of  central  disk  (nsc) 
if/  -  id  complex  yaw  and  pitch  of  central  disk  (nsc) 

y/  -id  complex  yaw  and  pitch  rate  of  central  disk  (nsc) 
n  =  3.4...nt  FEM  connectors  (nsc) 
n  =  3A...nt  FEM  connectors  (bfc) 
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5 


7rd2  /  4 

magnitude  of  projectile  velocity 
location  of  beam  ends 


V 


x0l,x23  dimensionless  length  of  fore  and  aft  aerodynamic  extensions 

xc  axial  location  of  center  of  mass 

a  angle  of  attack  of  central  disk  (nsc) 

ft  angle  of  sideslip  of  central  disk  (nsc) 


r 

Zu 

SE 

4> 

k 

K 

p 

pi 

a 

coi 

(Or 


— — ,  complex  cant  of  disk 

dx 

maximum  strain  of  rod 

SE  +  iSEz ,  lateral  displacement  of  disk  (nsc) 

roll  angle 

frequency  of  k-th  mode 
damping  of  k-th  mode 
air  density 

axial  variation  of  mass 
coJcqr 

lowest  elastic  frequency  of  beam  in  vacuum 
rigid  projectile  frequency  for  ,y,  =0,1 


F  =  ( Fx ,  Fy ,  F.  j  aerodynamic  force  exerted  on  missile  (nsc) 

M  =  ( M x ,  M y  ,M_  j  aerodynamic  moment  exerted  on  missile  (nsc) 

Fe{z}  real  part  of  z 

Im  {z}  imaginary  part  of  z 

Carat  superscript  denotes  quantity  for  a  single  element. 

Tilde  superscript  denotes  elastic  parameter  for  bent  missile. 
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B  subscript  denotes  parameter  for  bent  projectile. 

E  subscript  denotes  an  elastic  coordinate  parameter  (nsc). 
b  subscript  denotes  an  body-fixed  coordinate  parameter  (bfc). 


(bfc) 

body-fixed  coordinates 

(nsc) 

nonspinning  coordinates 
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